Model for frustrated spin-orbital chains: application to CaV 2 04 
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Motivated by recent interest in quasi-one-dimensional compound CaV2C>4, we study the ground 
states of a spin-orbital chain characterized by an Ising-like orbital Hamiltonian and frustrated in- 
teractions between S = 1 spins. The on-site spin-orbit interaction and Jahn- Teller effect compete 
with inter-site superexchange leading to a rich phase diagram in which an antiferro-orbital phase is 
separated from the orbital paramagnet by a continuous Ising transition. Two distinct spin liquids 
depending on the underlying orbital order are found in the limit of small spin-orbit coupling. In the 
opposite limit, the zigzag chain behaves as a spin-2 chain with Ising anisotropy. The implications 
for CaV2C>4 are discussed. 



Orbital degrees of freedom have been shown to play 
an important role in understanding the electronic and 
magnetic properties of transition metal oxidesi This is 
particularly so for frustrated magnets with partially filled 
orbitals. A well-studied case is vanadium spinels AV2O4 
where the A-site sublattice is occupied by divalent ions 
such as Zn 2+ or Mn 2+ , and the B-site magnetic V 3+ ions 
form a pyrochlore lattice.—^ It is known that geometrical 
frustration of classical spins on such a lattice precludes 
simple Neel ordering and gives rise to a highly degener- 
ate ground state. Orbital ordering in these compounds 
partially relieves the frustration by creating disparities 
in nearest-neighbor exchange constant, hence setting the 
stage for magnetic ordering at lower temperatures^^ 

In this paper, we investigate the physics of frustrated 
vanadium chains in which the interplay of geometrical 
frustration, spin-orbital couplings, Jahn-Teller effect, and 
enhanced quantum fluctuations leads to a rich phase di- 
agram. Our work is partly motivated by an attempt to 
understand another vanadium compound CaV204, which 
at room temperature crystalizes in the orthorhombic 
calcium- ferrite structure (space group Pnam) ,2iL2i2iI£ 
Contrary to its spinel cousins, V 3+ ions in CaV204 are 
arranged in zigzag chains of edge-sharing VOg octahedra 
(Fig. [TJ . Antiferromagnetic interaction on zigzag chains 
consisting of triangular loops is subject to geometrical 
frustration as well. The rather weak and frustrated inter- 
chain couplings make the vanadium chains quasi- ID sys- 
tems susceptible to quantum fluctuations. Couplings of 
vanadium orbitals to spins and phonons add yet another 
dimension to the intriguing physics of zigzag chains. 

In CaV204, the zigzag geometry results in a spin-1 
chain with comparable nearest and next nearest-neighbor 
interactions. Combined with the observation that 3d 2 
configuration of V 3+ ions tend to have an easy-plane 
anisotropy^ the quasi-lD compound CaV204 has been 
a favorable candidate for the long-sought chiral spin liq- 
uid where a long-range chiral order coexists with alge- 
braically decaying spin correlations ! 10 ! 11 However, as first 
pointed out by Pieper et al.^ orbital degeneracy in this 
compound drastically changes the above picture. Since 
the tig levels are split into a singlet and a doublet due to 
a flattened V06 octahedron, with the low-energy singlet 




FIG. 1: (Color online) (a) Two crystallographically inequiva- 
lent vanadium chains in CaV2 0zt. The V 3+ ions are arranged 
in zigzag chains of edge-sharing V06 octahedra. For both 
vanadium sites, the VO6 octahedra is flattened at room tem- 
peratures. A local reference frame is defined in such a way 
that z-axis is parallel to the tetragonal axis of the crystal field. 
The black and white circles denote the vanadium and oxygen 
ions, respectively, (b) A simplified view of the zigzag chain. 
The red, green, and blue bonds are parallel to the local [Oil], 
[101] and [110] axes, respectively. Consequently, electron hop- 
ping on the red, green, and blue bonds is possible only when 
dyz, d zx and d xy orbitals are occupied, respectively. 



always occupied, a double degeneracy remains for the 
other electron. The orbital degrees of freedom in V 3+ 
ion is thus described by an Ising-like variable. 

Recent experimental studies on CaV204 seem to rule 
out a possible chiral liquid phase as well. A structural 
transition at T s sa 141 K reduces the crystal symme- 
try to monoclinic P2i/n. 8,9 It was suggested in Ref. [1 
that the simultaneous orbital ordering relieves the mag- 
netic frustration of zigzag chains. As the inter-chain frus- 
tration is also lifted by the lattice distortion, a three- 
dimensional Neel order sets in at Tn sa 71 K. The 
collinear spins are parallel to the crystal 6-axis as evi- 
denced by both nuclear magnetic resonance and neutron 
diffraction measurements 

To make progress toward an understanding of the 
ground-state structure and the nature of phase transi- 
tions in CaV204, here we study the zero-temperature 
phase diagram of its building blocks, i.e. zigzag chains 
with S — 1 spins and Ising orbital variables. We pro- 
pose a theoretical model which includes the superex- 
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change (SE) interaction, relativistic spin-orbit (SO) cou- 
pling, and Jahn- Teller (JT) effect. We find that antiferro- 
orbital order favored by an Ising-like orbital exchange 
is destroyed in the presence of large on-site spin-orbit 
or Jahn- Teller coupling. Depending on the underlying 
orbital configuration, magnetic properties of the zigzag 
chain is equivalent either to those of two weakly coupled 
5 = 1 chains, or of an unfrustrated spin-1 ladder. In the 
limit of large spin-orbit coupling, the zigzag chain can 
be viewed as a spin-2 chain with anisotropic interaction. 
Finally, we discuss implications for CaV204. 

Model Hamiltonian. We first define a local reference 
frame for the two crystallographically inequivalent vana- 
dium chains (referred to as type-I and II chains here) in 
CaV2 04 such that the z-axis is parallel to the tetrag- 
onal direction of flattened VC>6 octahedron [Fig. Q^a)]. 
Nearest-neighbor bonds along the chain are parallel to 
the local [Oil] and [101] directions alternatively. In 
the following, we employ a convention in which even- 
numbered bonds are along the [Oil] axis. 

Geometrical frustration of the zigzag chain comes from 
the fact that second-nearest-neighbor bonds parallel to 
local [110] axis have a length close to that of nearest- 
neighbor bonds. In fact, the second-nearest-neighbor in- 
teraction is the dominant one as d xy orbital is occupied 
at every site due to the flattened VOe octahedra. The re- 
maining orbital degeneracy is described by a pseudospin- 
i with er z = ±1 corresponding to the \yz) and i\zx) 
states, respectively. 

Having introduced the basic notations, we now discuss 
a minimal model for the spin-orbital chain. Since orbital 
interaction with a 90° angle between vanadium-oxygen 
bonds is governed by direct dda exchange of t^g orbitals, 
the corresponding spin interaction on a bond depends on 
whether the relevant orbital is occupied^ Essentially, 
orbitals participate in the superexchange only via orbital 
projectors P yz/zx . Noting that P yz/zx = (1 ± cr 2 )/2, we 
first define the antiferro-orbital and ferro-orbital bond 
operators 

with ± signs for even and odd numbered bonds, re- 
spectively. The model Hamiltonian is divided into three 
parts: H = Hq + H± + iJ on _ s i te . The first Hq term rep- 
resents decoupled spin and orbital systems 

Hq = Ji ^™ ' Sn+2 — + K'On,n + l)- (1) 

n n 

Since Ji couples every second-nearest neighbors, the spin 
part can be viewed as two decoupled S = 1 Haldane 
chains, corresponding to the two blue lines in Fig. [ljb). 
The K and K' terms denote the energy gain of an 
antiferro-orbital and a ferro-orbital bond, respectively. 
In general, we have K > K' due to a finite on-site 
Hund's coupling Jjy, hence an antiferro-orbital order in 
the ground state of Hq. 



The H± term introduces interactions between the two 
spin-1 chains: 

H± = ^^(Ji(5„ in +i — J(O n! „+i) S„ • S„ + i. (2) 

n 

As discussed above interaction between spins depends 
on orbital occupations: the antiferromagnetic coupling 
Ji is nonzero only when d yz (d zx ) orbital is occupied at 
both sites of an even (odd) bond, whereas the strength of 
ferromagnetic J[ term depends on the expectation value 
of the antiferro-orbital bond operator C„ in+ i. 

Explicit expressions relating the exchange constants to 
microscopic parameters are obtained from the SE Hamil- 
tonian of vanadium spinels, 1 - where neighboring VOq oc- 
tahedra share the same edge as in zigzag chains consid- 
ered here. Assuming an exact octahedral site-symmetry, 
we find K = (1 + 2ry)^, J 2 = Ji = K' = (1 - 77)^, and 

,/{ = rjjj to lowest order in Hund's parameter ij = Jh/U. 
Here t is the hopping integral, and U is the on-site 
Coulomb repulsion. The parameters of the model can 
be estimated from known values of the same parame- 
ters in other vanadium compounds. Measurements on 
cubic vanadates yield Jh — 0.68 eV, U ~ 6 eV, and 
t ~ -0.35 eV^iii which gives 77 ~ 0.11 and t 2 /U ~ 
20.4 meV. The estimate of A is less certain, we find 
A ~ 13 - 25 meV^^^^ In the following, we shall mea- 
sure energy in units of t 2 /U. 

The remaining single-ion interactions are included in 
the Hamiltonian 

-ffon-sitc = -^X] ™ 5 ™ ~~ 5 ^2 a n- ( 3 ) 

n n 

The first term originates from relativistic SO coupling 
A(L • S). Since d xy orbital is always occupied, the x 
and y components of the orbital angular momentum are 
quenched; the remaining L z = —a x in the pseudospin 
representation. A similar situation has been studied in 
cubic vanadates^ The effect of the monoclinic structural 
transition at T s is modeled by the second term with 25 de- 
noting the level splitting due to the induced orthorhombic 
distortion of VOg octahedra. Note that the orthorhombic 
distortion is different on type-I and II chains, 8 '? 

Orbital orders. We first consider a simpler case of 
the model Hamiltonian by assuming the presence of a 
collinear Neel order on the zigzag chain. This is a plau- 
sible assumption as the SO term Xa^S^ breaks the spin 
SU(2) symmetry and, as will be discussed later, closes the 
energy gap of longitudinal magnons at large enough A, 
hence signaling a transition to the Neel state with S* par- 
allel to ±z. Even with this simplification, the competi- 
tion between inter-site exchange and various on-site inter- 
actions still poses a rather nontrivial problem. This study 
also sheds light on orbital orders in the ground state 
of CaV204, where spins develop a three-dimensional 
collinear antiferromagnetic order at Tjv ~ 71 K. 

Due to the strong second nearest-neighbor exchange 
J2, collinear orders on a zigzag chain consisting of re- 
peated + H spins have a quadrupoled unit cell. There 
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FIG. 2: (Color online) Spin and orbital orders on a zigzag 
vanadium chain. Neel orders (a) and (b) are related to each 
other by lattice translations. There are a total of four de- 
generate Neel states; the other two are related to states (a) 
and (b) by time reversal, (c) antiferro-orbital order consisting 
of staggered d yz and d zx orbitals. (d) and (e) correspond to 
ferro-orbital orders with real orbital d yz and complex orbitals 
(d yz ± id zx )/y/2, respectively. 
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FIG. 3: (a) Schematic phase diagram of orbitals in the pres- 
ence of a collinear Neel order on the zigzag chain, effectively 
described by Hamiltonian @ . AF and PM refer to antiferro- 
orbitally ordered and orbital paramagnetic phases, respec- 
tively, (b) Critical boundaries A c and 8 C measured in units 
of t 2 /U versus Hund's parameter r\. The critical distortion is 
determined analytically from S c — (2K — K' — Ji)/2, whereas 
A c is obtained from DMRG calculation. The various effec- 
tive parameters in Eq. Q are expressible in terms of model 
parameters Ji, J[, J2, K, and K' of the original SE Hamil- 
tonian, whose relations to r\ and t 2 /U can be found in the 
text. 



are a total of four degenerate Neel states related to each 
other by lattice translations and time reversal [Figs. Ufa) 
and (b)]. After applying a 7r-rotation about z-axis to 
pseudospins at < sites, we obtain an effective or- 
bital Hamiltonian: 

fforb = £[J + (-l) n *K<H ~ 5>,*S + h x a x n ). (4) 

n n 

This model is equivalent to an Ising chain with al- 
ternating nearest-neighbor couplings in a skewed mag- 
netic field. The effective exchange constants are J = 
(2K - K')/4 and K = v 2 {2J[ + Ji)/4. Here we have as- 
sumed (S„ • Sft+i) = T^ 2 on even and odd bonds, respec- 
tively. The longitudinal and transverse fields are given by 
h z = 5 + v 2 Ji/2 and h x = A, respectively. The parameter 
v < I characterizes the magnitude of the Neel order. Its 
value can only be determined with a proper treatment of 
the SO coupling term. For simplicity, we set v = 1 in 
the following calculation. Hamiltonian (fj| without the 
staggered exchange JC is one of the simplest models ex- 
hibiting nontrivial quantum critical points numerical 
calculation shows an order-disorder transition belonging 
to 2D Ising universality class. The staggered exchange 
K, involves higher-order spatial derivatives in the contin- 
uum limit and thus represents an irrelevant perturbation 
in the RG sense. 

A schematic phase diagram of model Q is shown 
in Fig. [3fa), where an antiferro-orbital phase is sepa- 
rated from the orbital paramagnet by an Ising transi- 
tion line. Along the <5-axis (A = 0), the antiferro-orbital 
phase coexists with ferro-orbitally ordered phase, shown 
in Figs. [He) an d (d) respectively, at the multicritical 
point S c = (2K — K' — Ji)/2. On the other hand, in the 
large A limit, the pseudospins are polarized by the trans- 
verse field h x such that spins S z = ±1 are accompanied 



by complex orbitals -^=(d yz ± id zx ), respectively, result- 
ing in a uniform orbital occupation n yz = n zx = 1/2 
at all sites [Fig. [21(e)] . Using the infinite-system DMRG 
method with periodic boundary condition, 1 - we obtain 
a critical A c at 5 = 0, taking into account the effect of 
staggered JC. The Ising transition line is bounded by 
the critical points (<5 C ,0) and (0, A c ). Fig. Efb) shows 6 C 
and A c as a function of parameter r\ = Jh/U. As ex- 
pected, the antiferro-orbital phase shrinks with decreas- 
ing Hund's coupling 3 a- In a full quantum treatment of 
the zigzag chain, we expect a similar critical line charac- 
terized by massless orbital excitations (Fig. [4|. 

Spin liquid phases. We now discuss the original Hamil- 
tonian (|T|)— (J3j) in small A limit without assuming the exis- 
tence of a magnetic order. It is important to note that the 
spin part of the Hamiltonian in this limit preserves a con- 
tinuous SU(2) symmetry, which can not be broken in one 
dimension. One thus expects stable spin-liquid phases 
whose properties depend critically on orbital configura- 
tions. Furthermore, the absence of a x term at A = 
indicates that the orbital part is described by a classi- 
cal Ising-like Hamiltonian. Consequently, the search of 
ground states reduces to first enumerate over all possible 
Ising configurations {c^}, and then compare their en- 
ergies taking into account contribution from spins. As 
antiferro-orbital and ferro-orbital orders are favored by 
SE and JT interactions, respectively, it is natural to con- 
sider these two configurations first. 

In the case of antiferro-orbital order [Fig. [He)] where 
bond operators (O n , n +i) = and (C„ !?1 +i) = 1, the 
zigzag chain behaves as two spin-I chains weakly coupled 
by a frustrated ferromagnetic J[ . The corresponding spin 
liquid phase (SL1 phase in Fig. [4|) has an energy gap at 
an incommensurate wavevector On the other hand, the 
frustrated J[ coupling is quenched by the ferro-orbital 
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FIG. 4: (Color online) Schematic phase diagram of the zigzag 
vanadium chain. The various phases are characterized by the 
magnetic properties. SL and N represent spin liquid phase 
and Neel order, respectively. The SL1 and Nl phases are 
accompanied by an antiferro-orbital order, while in the N2 
and SL2 phases the orbitals behave as an orbital paramagnet. 
The first-order transition along the dashed line is an extension 
of the multicritical point (8 C , 0) in the phase diagram Fig. 



order shown in Fig. Eld). Depending on the sign of ferro- 
Ising order (ofj = ±1, the Jj term is nonzero only on 
even or odd bonds, but not both. The spin Hamiltonian 
is equivalent to that of a spin-1 ladder with rung cou- 
pling J\. The magnetic ground state is again disordered 
(SL2 phase in Fig. H])i 21 i 22 The magnetic energy of spin-1 
ladder with arbitrary rung coupling has been calculated 
using quantum Monte Carlo method in Rcf. Com- 
paring the energy of the two phases, including both spin 
and orbital contributions, yields a boundary S c surpris- 
ingly close to the one obtained assuming a preexisting 
Neel order. 

Although the SU(2) symmetry is broken in the pres- 
ence of SO coupling, both spin liquids persist up to fi- 
nite A. For SL2 phase stable at 6 > S c , SO coupling 
provides an easy-axis anisotropy DS% with D ~ —X 2 /26. 
The spin-1 ladder undergoes an Ising transition to a Neel 
phase (N2 phase in Fig. [4]) with increasing A. At large S, 
the condition D = D c gives rise to a critical A c oc vo. On 
the other hand, Neel transition at small 8 can be under- 
stood by considering the A dependence of the spin gap in 
SL1 phase. The singlet ground state of weakly coupled 
Haldane chains has triply degenerate magnon excitations 
with dispersion uik — + v 2 (k — fco) 2 - In the limit of 

vanishing J[, the system reduces to two decoupled Hal- 
dane chains with fco = ir, Ao = 0.410J2 (Haldane gap for 
spin-1), and v = 2.49^22- A simple second-order pertur- 
bation calculation adds a correction to the spin gap 
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where A CT m 2(5 C — S) is the energy of a domain- wall 
pair (a flipped pseudospin in the antiferro-orbital state), 
|0) is the singlet ground state, and \k) denotes one- 
magnon excitation with wavevector fc. The matrix el- 
ement (0\S*\k) = \fZe ikn , with Z w l. 23 For A CT > A , 
the spin gap decreases with increasing A and eventually 



reaches zero at A C 2 oc y/2(5 c — S) — Ao, signaling a tran- 
sition into magnetically ordered phase (Nl phase). 

The two Neel phases in Fig. 2] are distinguished by 
the underlying orbital orders, the phase boundary A c i 
separating Nl from N2 phases is thus analogous to the 
Ising transition line of the orbital only model [Fig. EJb)]. 
Since the upper critical A c i oc Jh, at small Hund's cou- 
pling the zigzag chain could bypass the Nl phase and en- 
ter the orbital paramagnetic phase simultaneously with a 
magnetic ordering. A conjectured phase diagram of the 
zigzag chain is shown in Fig. 0J 

Anisotropic J — 2 spin chain. In the limit of large SO 
coupling, the appropriate degrees of freedom are effective 
spins of length J = 2, where J = L' + S and L' = 1 is 
the angular momentum of the unoccupied ti g hole. The 
system thus behaves as a spin-2 chain with anisotropic 
exchange interaction^- Furthermore, JT coupling adds an 
anisotropy DJ 2 + E(J 2 — J 2 ), where D < and E oc 
5 are proportional to the tetragonal and orthorhonibic 
distortions of V06 octahedron, respectively. Assuming a 
dominating easy-axis anisotropy \D\ ^> E, the N2 phase 
can also be viewed as Neel ordering of the effective spins 
J z = ±2 in such a way that spins of a given direction 
S z = ±1 is coupled to orbital angular momenta L' z = ±1, 
respectively. The corresponding ferro-orbital order with 
complex orbitals d yz ± id zx is shown in Fig. [2je). The 
V 3+ ions have a reduced magnetic moment fi = (2S — 
L')fi B = lpL B . 

Discussion. We have presented and analyzed a min- 
imal model of frustrated vanadium chains, pertinent to 
quansi-lD compound CaV204. A conjectured phase dia- 
gram (Fig. 3]) is obtained based on analytical arguments 
and numerical calculations. The observed P2±/n crys- 
tal symmetry of CaV204 at low temperatures indicates 
that only two inequivalent vanadium sites exist as in the 
high temperature phase. The absence of doubled unit 
cell resulting from the antiferro-orbital order thus implies 
that both vanadium chains are in the orbital paramag- 
netic phase. This is a plausible conclusion noting that a 
rather small A c i w 0.22 (t 2 /U) m 4.5 meV is estimated 
from Fig. [3] assuming rj ss 0.11. However, the monoclinic 
distortion below T s places the two types of vanadium 
chain at rather distinct regions of the phase diagram. 

The type-I chain acquires a small 5 in addition to 
the dominating tetragonal crystal field and behaves as 
a spin-2 chain subject to an easy-axis anisotropy. This 
Ising anisotropy is important to the stabilization of three- 
dimensional Neel order in CaV204, as the collinear spins 
are found to point along the easy axis of type-I chains. 8 
The measured V moment 1.06/is is also consistent with 
the picture of an anisotropic spin-2 chain>£ On the other 
hand, a strong orthorhombic distortion S at type-II ions 
completely removes the orbital degeneracy and makes the 
vanadium chains effectively spin-1 ladders. In fact, the 
well separated t 2g levels at type-II ions lead to a pos- 
sible easy-plane spin anisotropy. 8 Consequently, induced 
collinear order on type-II chains follows the spin direction 
of type-I vanadium ions, as indeed observed in CaV204. 
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